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Systematic analysis of the deformations of proton and neutron densities in even-even C isotopes 
was done based on the method of antisymmetrized molecular dynamics. The E2 transition strength 
was discussed in relation to the deformation. We analyze the B{E2\ 2^ — > 0^) in ^®C, which has been 
recently measured to be abnormally small. The results suggest the difference of the deformations 
between proton and neutron densities in the neutron-rich C isotopes. It was found that stable proton 
structure in C isotopes plays an important role in the enhancement the neutron skin structure as 
well as in the systematics of B{E2) in the neutron-rich C. 



I. INTRODUCTION 



In light unstable nuclei, exotic phenomena such as neutron halo and neutron skin structures were discovered owing 
to the progress of the experimental technique. These contradict to the traditional understanding for stable nuclei 
where the proton density and neutron densities are consistent with each other in a nucleus. These phenomena imply 
On ] that the exotic features may appear in unstable nuclei due to the difference between proton and neutron densities. 
, Another subject concerning the difference is the deformations of proton and neutron densities. For example, the 
opposite deformations between proton and neutron densities in proton-rich C isotopes were theoretically suggested 

^ : [1]. 

• Recently, the life time of the 2f state of ^^C has been measured [2]. It indicates the abnormally small E2 transition 
; strength as B{E2; 2+ ^ 0+) =0.63 e^fm" in ^^C, compared with those for other C isotopes(i°C, ^^C and ^''C). As well 
■ known, B{E2) is related to the intrinsic deformation of the nucleus. Considering the excitation energy Ex{2f) = 1.766 
00 ' MeV of ^®C, it is expected that this nucleus is not spherical but has a deformed structure. In case of normal stable 
. nuclei, B[E2) values are generally large in the deformed nuclei. It means that the hindrance of the B{E2; 2^ 0^) in 
' ^^C seems to contradict to the deformation expected from the Ej;{2^). Moreover, the neutron and proton transition 
,. matrix elements, Mn, Mp derived from the ^°*Pb-|-^^C inelastic scattering [3] imply that the neutron excitation is 
' dominant in the 2^ state of ^^C. From the point of view of collective deformations, these experimental results suggest 
I ] the possible difference between proton and neutron shapes in ^^C. 

Q ■ In the theoretical work on the proton-rich C isotopes [1] , it was suggested that the difference of the shapes between 
^ ' proton and neutron densities may cause the suppression of B{E2). In the previous work [1], the difference between 
' the proton and neutron shapes in the neutron-rich C isotopes was predicted as well as the proton-rich side. It is 
' natural to consider that the systematic analysis of the proton and neutron deformations in connection to the B{E2) 
• ^ is important to understand the properties of the neutron-rich C. It may be helpful also to predict exotic phenomena 
in unstable nuclei. 

^ , In this paper, we study the deformations of proton and neutron densities in C isotopes based on the theoretical 
- - ' calculations with antisymmetrized molecular dynamics (AMD). The AMD method is a useful approach for the structure 
study of stable and unstable nuclei. The applicability of this method in the systematic analysis of the light nuclei 
has been proved in many works [4-8] . Especially, in the description of the deformation and clustering aspect in light 
nuclei, this method has an advantage over the normal mean-field approaches like Hartree-Fock calculations, which 
sometimes fail to describe the deformation in the very light nuclei such as ^^C. Although shell model calculations are 
useful to investigate the level structure of light nuclei, they are not suitable for the systematic study of the deformations 
because of the difficulty in directly extracting intrinsic deformations as well as the ambiguity of the effective charges 
in the unstable nuclei. In order to extract a naive picture on the proton and neutron shapes we apply the simplest 
version of AMD, based on a single AMD wave function. By using the AMD method, we analyze the systematics of 
the deformation and E2 transition strength in neutron-rich C isotopes. The hindrance of the E2 transition strength 
in ^^C is discussed. The theoretical predictions for the deformation and B{E2) in further neutron-rich isotopes, ^^C 
and ^°C, are also reported. 

This paper is organized as follows. In section II, the formulation of AMD is briefly explained. We show the results of 
the energies, radii, and B{E2) obtained by the simple version of AMD, while comparing them with the experimental 
data in Sec. IV. In Sec.V, the intrinsic deformations of proton and neutron densities are analyzed in connection to the 
observable such as the E2 and radii. In Sec. VI, we show the results of ^^C obtained by an extended version of AMD. 
Finally, a summary is given in Sec. VII. 
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II. FORMULATION 



The formulation of AMD for nuclear structure studies is explained in [4,6,9]. 

The wave function of a system with the mass number A is written by a superposition of AMD wave functions 
^AMD- An AMD wave function is given by a single Slator determinant of Gaussian wave packets as, 

$AMd(Z) = -^A{ip-L, ip2, </?a}, (1) 

where the i-th single-particle wave function is written as, 

= 4':x-^XiT^, (2) 
<^Xi(rj) oc exp{-u{rj - ^f}, (3) 

Xi = i^+^i)Xf + i^-^i)Xl- (4) 

In the AMD wave function, the spatial part is represented by complex variational parameters, Xii, X2i, Xa^, which 
indicate the center of the Gaussian wave packets. The orientation of the intrinsic spin is expressed by a variational 
complex parameter ^j, and the iso-spin function is fixed to be up(proton) or down(neutron). 

In order to obtain a naive understanding on the systematics of intrinsic deformations, we use a simplest version of 
the AMD method which was applied for Li, Be and B isotopes in Ref. [4]. Namely, we perform energy variation for a 
parity-eigen state, P'^^amd = ^amd^ projected from an AMD wave function. We consider the AMD wave function 
obtained by the energy variation as the intrinsic state, and the total-angular- momentum projection(P^^) is done 
after the variation to evaluate observables such as the energies, radii, and the transition strength. Thus the variation 
is done after the parity projection, but the total-angular-momentum projection is performed after the variation. This 
method is called as VBP(variation before projection) in the present paper. For further investigations of the level 
scheme of ^^'C, we also perform the VAP(variation after projection) calculation with respect to both the parity and 
total-angular-momentum projection as the same way as done in Refs. [9,10]. In the VBP calculations, we fix the 
orientation of the intrinsic spin to be up or down. In the VAP calculations, ^j's are treated as free variational 
parameters. 



III. INTERACTIONS 



The effective nuclear interactions adopted in the present work consist of the central force, the spin-orbit force and 
the Coulomb force. We adopt MVl force [11] as the central force. This central force contains a zero-range three-body 
force as a density-dependent term in addition to the two-body interaction. The Bertlett and Heisenberg terms are 
chosen to be & = /i = 0. We use the parameter set, case 3 of MVl force with the Majorana parameter as m = 0.576, 
which was adopted in Ref. [4]. Concerning the spin-orbit force, the same form of the two-range Gaussian as the GSRS 
force [12] is adopted. The strengths of the spin-orbit force, (a) uj = —uu = uis = 900, and (b) 1500 MeV are used. 
In the VAP calculations, we also use the interaction parameters, (c) the case 1 of MVl force with m = 0.62 and the 
spin-orbit force with uis = 3000 MeV used in Ref. [9]. 



IV. RESULTS OF VBP CALCULATIONS 



The structure of positive parity states of even-even C isotopes are studied by the VBP calculations within the 
framework of AMD. In this section, we present theoretical results, such as the energies, radii, and E2 transitions, 
while comparing them with the experimental data. The optimum width parameter v in Eq.3 is chosen to mini- 
mize the energy of the system for each nucleus. The adopted v parameters are listed in Table. I. After the varia- 
tion, we perform the total-angular-momentum projection Pjj^^ and diagonalize the Hamiltonian and norm matrices, 
(^'mk'*amdI^I^mk"*amd) and (P^'^x'*AAfDl^'MK"*^MD) with respect to the i^-quantum {K',K"). Conse- 
quencely, we obtain 0^ and 2f states which can be regarded to belong to the ground K'^ = 0+ band in each C 
isotope. In cases of ^°C, ^^C and ^^C, the second 2+ state in the K'^ = 2+ band arises because of the axial asymme- 
try. 

The binding energies and the excitation energies of the low- lying positive-parity states are shown in Fig.l and in 
Table.II. The binding energies of C isotopes are systematically reproduced by the present calculations (Fig.l). As 
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TABLE I. The a(loj)tcd width parameter ;/ of the AMD wa\'e functions for C isotopes. 
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FIG. 1. Binding energies of C isotopes. Solid(dotted) hne indicates the VBP calculations with m = 0.576 and uis = 1500(900) 
MeV. 



mentioned before, the second 2+ state, which belong to the side band = 2+, is found in addition to the 2]*" state, 
in IOC, and ^^C. 

The excitation energies of 2'^ states tend to be underestimated by the calculations, especially in the neutron- 
rich C isotopes. This is a general tendency of the VBP calculations where the 0]*" and 2^ states are obtained by the 
total-angular-momentum projection from a single intrinsic wave function. It is considered to be because the wave 
fimction obtained in the VBP may be optimized for the 2^ state while it is relatively not enough for the 0^ state to 
reproduce the level spacing between the 2^ and 0^ states. In other words, the energy of 0^ state is considered to 
be relatively overestimated in the VBP. The quantitative reproduction of Ex{2'\^) is improved by VAP calculations, 
where the wave function is optimized for each spin-parity state. Based on the VBP and VAP calculations, the level 
structure of ^^C is discussed later in Sec. VI. We give a comment that also the AMD+GCM (generator coordinate 
method) calculations,which is an extended version of AMD, well reproduce the £'x(2]'") in neutron-rich C. 

Figure 2 shows the results of the root-mean-square radii comparing with the experimental radii derived from the 
interaction cross sections. In the systematics of the matter radii in the even-even C isotopes, there is a gap between 
^^C and ^^C. This behavior is reproduced by the present calculations and can be described by the large deformation 
in ^^C and the spherical shape due to the shell closure in ^^C. The details are discussed in the next section. 



TABLE II. Excitation energies of 2+ states of C isotopes obtained by the VBP calculations. 
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FIG. 2. Root-mean-square radii of C isotopes. The matter radii are displayed in the upper panel(a). The solid(dotted) line 
indicates the VBP calculations with m = 0.576 and = 1500(900) MeV. The experimental data, which are derived from 
the interaction cross sections [13], are shown by cross points. The point-like proton and neutron radii calculated with with 
m = 0.576 and uu = 1500 MeV are shown in the lower panel(b). 
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TABLE III. E2 transition strength in C isotopes calculated by the VBP. The unit is e'^fin^. The experimental data for ^^C 
is taken from Ref. [2] and those for the other C are from Ref. [14]. 
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The results for the E2 transition strength are listed in Table. III. In the calculations, it is found that B{E2; 2^ Of) 
drastically changes with the increase of neutron number in C isotopes. The theoretical E2 strength in ^^C is very 
small, which is consistent with the recent measurement of the abnormally small B{E2) in ^^C [2]. The present results 
predict that, also in ^^C, B{E2; 2^ ) is very small as the same order of that in ^^C. It is interesting that 

the calculated E2 transition strength in ^"^C is not as small as those in ^^C and ^*C. The systematic change of the 
B{E2-,2^ Q'l) can be understood by the deformations of proton and neutron densities in the intrinsic states as 
described in the next section. 



V. DISCUSSION 



In this section, we analyze the intrinsic deformations of proton and neutron densities and discuss their effect on the 
observables such as the E2 transitions and radii. 



A. intrinsic deformation 



We display in Fig. 3 the deformation parameters (/?, 7) for the proton and neutron densities, which are defined by 
the moments {x^)^ {y'^)^ ^ind i^"^) of the intrinsic AMD wave function as. 



((x2>(j/2)(^2))l/6 




(5) 
(6) 
(7) 



Here the x, y, and z directions arc chosen so as to satisfy (.t^) < {]/) < (z^) and (xy) = (yz) = (zx) = 0. As seen in 
Fig. 3, we find the drastic change of the neutron deformation in C isotopes with the increase of the neutron number. 
The neutron deformations are prolate, oblate and spherical in ^"^C, ^^C and ^''C, respectively. In the neutron-rich 
region, it become prolate again in ^*'C, and they are triaxial and oblate in and ^°C, respectively. In contrast to 
the variation of neutron deformation, the proton deformations are rather stable. The deformation parameters for the 
proton densities lie in the oblate region 7^1- These behavior are the same as the results of Ref. [1]. 

By analysing the component of the _fi'-quantum states(P/^j^-$^^^jj) in the 2+, it is found that the 2+ stat c can be 
approximately written by a single K = state when we chose a proper axis. Then, we can define the (approximate) 
principal axis Z in the body-fixed frame and form the gound K = band with the and 2^1 states in each C isotope. 
In the systems, ^^C, ^''C and ^°C, with the oblate or spherical neutron shapes, the principal axis Z is the same as the 
symmetric axis x which has the smallest moment (x^) as shown in Fig. 4(b). In other words, the dominant component 
of the excited state 2^ are the Jz = K — state with respect to the symmetric axis x. It is consistent with a naive 
expectation for the collective rotation. It is notable that, in ^"C, ^^C and ^^C, the deformations are different between 
proton and neutron densities. In these nuclei, the symmetric axis for the proton shape differs from that for neutron 
density. Namely, the symmetric axis of the oblate proton density is the x direction, while that of the prolate neutron 
is the z direction. The schematic figure of the proton and neutron shapes in ^^C is illustrated in Fig. 4 (a). Such the 
configuration of the proton and neutron shapes is energetically favored because it has the maximum overlap between 
the proton and neutron densities. Because of the coexistence of the different shapes between proton and neutron, the 
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FIG. 3. Deformation parameters /3, 7 of the intrinsic states calculated with m = 0.576, uu = 900 MeV (a,b) and 1500 MeV 
(c,d). The mass numbers A are written by the corresponding points. 



second 2'^ state appear due to the triaxiality of the total system. By analysing the component of the K-qaaxitmn 
states {Pmk^amd) the excited states, it is found that the Of and 2f states belong to the ground = 0+ band 
and the state is classified into the side-band = 2+ when we regard the z direction with the largest moment 
(z^) as the principal axis Z as noted in Fig. 4(a). It is important that the principal axis Z is not the same as the 
symmetric axis x for the proton density but is perpendicular to the x in these nuclei. The existence of the side band 
= 2+ in ^^C and ^^C has not been experimentally confirmed yet. Concerning ^"C, the triaxiality of the mirror 
nucleus ^°Be were discussed in Ref. [15], and the known 23 state in ^°Be is assigned to be the band-head state of the 
side band K"" = 2+ [15,10]. 

We compare the matter deformation of the present results with a HF+BCS calculation with Skyrme force by Tajima 
et al. [16]. In the HF+BCS calculation, the C isotopes in the A < 14 region have spherical shapes, which contradict 
to the knowledge that ^^C is oblately deformed. It is natural because the mean-field calculation is considered not 
to be valid for very light nuclei. In the neutron-rich C, the calculated quadrupole deformation parameter is positive 
in ^^C and ^*C and oblate in ^°C according to the HF+BCF calculation. The general behavior of the quadrupole 
deformation of the matter density in neutron-rich C seems to be similar with the present results, though the different 
shapes between proton and neutron hardly appear in the HF calculations. 



B. E2 transition 



The intrinsic deformation is closely related to the E2 transition strength. Since the present results indicate that the 
proton radius does not drastically change with the increase of the neutron number in C isotopes as shown in Fig. 2, 
the B{E2) is dominantly determined by the deformations. 

As mentioned before, the Of and 2^ states belong to the ground = 0+ band, where we regard the x axis as the 
principal axis Z in ^^C, ^'^C and ^°C, and the z axis as the Z axis in ^°C, ^^C and ^*C. In order to link the intrinsic 
deformations with the B{E2), we remind the reader the well-known approximate relation between the B{E2) and the 
intrinsic quadrupole moment Qo- 

B{E2;2t ^Ot) = ^e^Ql (8) 
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FIG. 4. Schematic figures for intrinsic deformations of tlie proton and neutron densities, (a) tfie oblate proton and prolate 
neutron shapes in ^''C, and (b) the oblate proton and oblate neutron densities in ^"C. The x, y, and z axes are chosen to be 
{x^} < (j/^> < {2^>- The principal axis Z for the ground band Jz = K"" = is also displayed. 



The intrinsic quadrupole moment Qo here is defined with respect to the principal axis as Qo = 2(Z^) — (X^) — (F^) 
, and is related to the deformation parameter 7p for the proton density. In ^^C, ^^C and ^°C, where Z pa x and 
7 « 7r/3, Qo is approximated as: 

Qo = 2(,r2) - {y') - iz') « -^N,eP,rl (9) 

where Np and Ve are the proton number and the root-mean-square charge radius. 
In ^°C, ^^C and ^^C with Z z, Qo depends on the deformation parameter 7^ as: 

Qo = 2{z^) - {x^) - « ^J^Npepp cos^y,. (10) 

The important point is that the effect of the proton deformation on the Qo decreases in these nuclei because of 
the factor cos7p. Especially, in -'^^C and ^®C, the deformation parameter 7^ ~ 7r/3 makes the Qo very small. This 
is the reason for the unusually small B{E2; 2^ ) in ^^C and ^^C shown in Table. III. In other words, the 

B{E2; 2^ — > 0^) in ^^C and ^^C is reduced due to the deviation the principal axis Z of the total system from the 
symmetric axis x for the proton density. The origin of the deviation is the different shapes between proton and 
neutron densities. On the other hand, it is interesting that the larger B{E2 : 2^ Of) is predicted in ^°C, because 
it has the oblate proton and neutron shapes, therefore, the principal axis aligns to the symmetric axis x as shown in 
Fig.4(b). 

The hindrance of the B{E2 : 2^ ^ 0^) in ^^C can be described also from the point of view of the collective rotation. 
The 2^ state of ^^C is roughly regarded to be the state with = K = Q. The J = 2, = state is given by the 
linear combination of the = 2 state and the Jy = 2 state, which are rotating around the x and y axes, respectively. 
The rotation around x axis causes no proton excitation, therefore, gives no contribution on the B(E2), because x is 
the symmetric axis of proton density. As a result, the J^, = 2 component reduces the B{E2 : 2^ ^ 0^). In contrast 
to the small B{E2 : 2^ ^ 0^), it is predicted that the B{E2 : 2^ — » ) in the side band is large because the 23 
state is dominated by the Jz =2 state, where the oblate shape of the proton densities gives proton excitations. 

As mentioned above, the abnormal small B{E2 : 2^ Of) in -'^^C can be understood by the difference between 
oblate proton and prolate neutron shapes. Even if the proton shape in ^^C is spherical or slightly prolate, the small 
B{E2 : 2f 0^) can be also described. The characteristic of the present result is the stable proton structure in 
the series of C isotopes and the prediction of the second 2+ state with large B{E2) in the side band. In order to 
conclude the intrinsic shapes of proton and neutron densities, we need further experimental information such as the 
systematics of B{E2) in other neutron-rich nuclei and some probes for the side-band in ^^C. 



C. radii 



In Fig. 2, we present the calculated results of matter, proton and neutron radii of C isotopes while comparing them 
with the experimental radii which are derived from the interaction cross sections [13]. As seen in Fig. 2(b), the proton 
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FIG. 5. Level scheme of the low-lying states of C. The theoretical results obtained by the VAP and VBP calculations with 
the interaction parameters (b)case 3 of MVl with m = 0.576, uu = 1500 MeV and (c)case 1 of MVl with m = 0.62, uis = 3000 
MeV are illustrated with the experimental data. 

radius does not drastically change with the increase of neutron number, while the neutron radius increases rapidly in 
the neutron-rich ^ > 16 region. The gap of the matter radii between ^''C and ^^C is fomid to originate from a gap 
in neutron radii. The reason for the gap is described by the neutron deformations as follows. The neutron density 
is compact in ^^C because it has a spherical shape due to the neutron shell closure. On the other hand, in ^^C, the 
neutron radii is large because of the prolate deformation. 

The calculated matter radii systematically reproduce the experimental data. The large matter radii in the neutron- 
rich region originate from the enhancement of the neutron radii. In contrast to the variation of the neutron radii, the 
proton radii is stable and is generally compact. As a result of the stable proton structure, the neutron skin structure 
enhances in the neutron-rich C isotopes. The development of neutron-skin in the C isotopes found in the present 
results is consistent with the results of mean-field calculations [16]. In the calculations by Thiamova et al. [8] with 
AMD+GCM, the radii of neutron-rich C arc somehow underestimated. This is considered to be because of the lack 
of the three-body term in the effective force in their calculations. 

VI. RESULTS OF VAP CALCULATION 

So far, we discuss the structure of C isotopes based on the VBP calculations. As mentioned before, the quantitative 
reproduction of the excitation energy Ex{2^) in the VBP calculations is not satisfactory in the neutron-rich C. Instead, 
the VAP calculation [9,10] is more useful to describe the detail of the level spacing. We perform energy variation after 
the spin-parity projection for the , 2^ and 22 states of ^^C, and obtain three independent AMD wave functions. 
We evaluate the observables by diagonalizing the Hamiltonian and norm matrices with respect to the three wave 
functions as done in Refs. [9,10]. 

Figure 5 shows the level scheme of the low- lying states of ^^C obtained by the VBP and VAP calculations with 
the interaction parameters (b) case 3 of MVl with m = 0.576, uis = 1500 MeV. We also show the VAP results 
obtained by the other parameter set (c) case 1 of MVl with m = 0.62, uis = 3000 MeV, which were adopted in Ref. 
[9] to reproduce well the level structure of ^^C. The level spacing between O]*" and 2^ is well reproduced by the VAP 
calculations. The excitation energies of the side-band = 2+ states rise comparing those with VBP. Especially, in 
the VAP with the interaction (c), the 22 state becomes relatively high because of the strong spin-orbit force. 

The results of the E2 transition strength in ^^C are listed in Table. IV. Comparing the theoretical values with the 
experimental data, the VAP calculations tend to overestimate B{E2;2^ 0^). It should be noted that the triaxial 
shape 7p ~ 6/7r of the proton density is found in the VAP results. This is the reason for the larger B{E2; 2^ Of ) 
in VAP than that in VBP in which ^^C has an oblate proton shape 7p ~ 3/7r. Because of the triaxial proton shape 
in VAP, the E2 transition in the side band, B{E2; 2^ Of), is smaller than that of VBP. 

In the present VAP calculations of ^^C, the Of and 2^ level spacing is well reproduced, while the reproduction of the 
small B{E2; 2^ — > O^") is not satisfactory. In order to insight the structure of ^^C, we need to improve wave functions 
in more detail. Moreover, the further experimental information for other excited states is helpful to determine the 
proton shape in ^^C. 



TABLE IV. The E2 transition strength in C. The theoretical results are obtained by the VAP and VBP calculations with 
the interaction parameters (b) m = 0.576, uu = ISOOMeV and (c) m = 0.62, uu = SOOOMeV. The unit is e^&n*. 
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VII. SUMMARY 

We studied the structure of even-even C isotopes with the AMD method. The experimental data of the binding 

energies, B{E2), and radii of C isotopes arc well reproduced in the present calculations. Systematic analysis of the 
proton and neutron shapes in C isotopes was done based on the VBP calculations of AMD. The results indicate that 
the neutron shape drastically changes with the increase of neutron number, while the proton shape is rather stable. 
It is suggested that the difference between proton and neutron shapes may appear in ^^C and -'^^C as well as ^"C. The 
E2 transition strength, B{E2), was discussed in relation to the deformation. The unusually small B{E2; 2^'" ^ Oj'") in 
^^C, which has been recently measured, was described by the coexistence of the oblate proton shape and the prolate 
neutron shape. According to the present prediction, the B{E2; 2^ 0^) in ^^C is as small as that in -'^^C, while the 
B(E2) is larger in ^''C. The deviation between the proton and neutron shapes play an important role in the small 
B{E2). The present results show the enhancement of the neutron skin structure in neutron-rich C. It was found 
that the stable proton structure in C isotopes plays an important role in the neutron skin structure as well as in the 
systematics of B{E2). 

In order to extract a naive picture on the proton and neutron shapes we applied the simplest version of AMD, based 
on a single AMD wave function. The quantitative reproduction of both the excitation energy (2^)) and B{E2) was 
not satisfactory in the present work. We consider that more detailed investigations with the improved wave functions 
as well as with the use of appropriate effective nuclear forces are required. 

In the experiment of the -^^^Pb+^^C inelastic scattering, the contributions of nuclear excitation and Coulomb ex- 
citation from the ground state to the 2+ (1.77 MeV) state were analyzed [3]. The ratio of the neutron and proton 
transition matrix elements M„/Mp implies the neutron excitation is dominant in the 2f state in -"^^C. This is con- 
sistent with the present results of oblate proton and prolate neutron shapes. Unfortunately, with this experimental 
information for the excitation to the 2^(1.77 MeV), it is difhcult to know whether the proton deformation is oblate or 
slightly prolate (or slightly triaxial), because in both cases the sign of the Mp(Of ^ 2f ) is same. Namely, even if the 
proton density has an oblate shape, the Mp is not negative but a positive value because of the different orientations 
of the symmetric axes between the proton and neutron shapes as shown in Qo given by Eq.lO. The characteristic of 
the present result is the stable proton structure and systematics of B{E2) in the series of C isotopes. If the ^^C has 
the oblate proton shape, the second 2+ state appears. We should stress that, whether the proton density is oblate or 
slightly prolate in -'^^C, the proton shape must differ from the large prolate deformation of the neutron to decribe the 
small B{E2). It is concluded that the small B{E2) indicates the difference between proton and neutron shapes in ^®C. 
In order to understand the details of the intrinsic shapes of proton and neutron densities, we need more systematic 
analysis of C isotopes with the help of the further experimental information such as the B{E2) in other neutron-rich 
nuclei, and ^"C, and information for the side-band in ^®C. 
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